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Abstract 

Directed and elliptic flow for the 197 Au + 197 Au system at incident energies 
between 40 and 150 MeV per nucleon has been measured using the INDRA 4vr multi- 
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detector. For semi-central collisions, the elliptic flow of Z < 2 particles switches 
from in-plane to out-of-plane enhancement at around 100 MeV per nucleon, in 
good agreement with the result reported by the FOPI Collaboration. The directed 
flow changes sign at a bombarding energy between 50 and 60 MeV per nucleon and 
remains negative at lower energies. The conditions for the appearance and possible 
origins of negative flow are discussed. 

Key words: Symmetric heavy ion collision, Squeeze-out, Directed flow, Elliptic 
flow 
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Considerable progress has been made recently in determining the equation 
of state of nuclear matter from heavy-ion reaction data [1,2]. A prominent 
role among the available observables is played by the collective flow as it is 
most directly connected to the dynamical evolution of the reaction system, 
including the momentum dependence of nuclear interactions and in-medium 
effects (for reviews see [3,4]). Very significant constraints on the possible range 
of interaction parameters have been derived from transverse and elliptic flow 
variables [2]. 

For 197 Au + 197 Au collisions, the amplitudes of transverse and elliptic collec- 
tive motion assume their maxima at bombarding energies of 300 to 400 MeV 
per nucleon. At these energies, the sign of the elliptic flow indicates a pref- 
erence for emissions perpendicular to the reaction plane. The change-of-sign 
recently observed at ultra-relativistic energies has received particular interest 
as it reflects the increasing pressure buildup in the non-isotropic collision zone 
[5,6]. Toward the lower incident energies, the directed flow observables are 
presumed to be related to the competition of mean-field and nucleon-nucleon 
collision dynamics and their evolution with the bombarding energy [7,8,9]. 
Elliptic flow has been identified with the collective motion resulting from the 
rotation of the compound system or the expansion of the hot and compressed 
participant zone, possibly modified by the shadowing effect of the colder spec- 
tator matter [10,11,12,13,14,15]. Also here, at the intermediate energies, the 
transition energies at which the flow parameters change sign are particularly 
useful for the comparison with theory. Their correct prediction requires the 
cancellation of the competing momentum components which is highly sen- 
sitive to specific parameters of the theoretical description. A more technical 
advantage is the weak sensitivity of the transition energies to the resolution 
achieved in reconstructing the reaction plane. 

In this Letter, we present results of the flow analysis applied to the data for 
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197 Au + 197 Au collisions at incident energies from 40 to 150 MeV per nucleon, 
collected with the 4-rr INDRA multi-detector [16] and with beams provided 
by the heavy-ion synchrotron SIS at GSI. With this energy range, the gap 
is bridged between existing excitation functions of collective motion in the 
relativistic and Fermi-energy domains. In agreement with the data of other 
groups, the transition from predominantly in-plane to out-of-plane emissions 
at about 100 MeV per nucleon, as reported by the FOPI collaboration [15], 
is confirmed. Directed flow is found to change its sign at a bombarding en- 
ergy below 60 MeV per nucleon [8], however with parameters that are found 
to depend strongly on the exact method applied and on the experimental 
acceptance. 

Details of the experiment, including the identification and calibration proce- 
dures have been presented in Refs. [17,18] and references given therein. For 
impact parameter selection, the total transverse energy E± of light charged 
particles (Z < 2) was used as a sorting variable. The minimum-bias distribu- 
tions of this quantity scale in proportion to the center-of-mass collision energy, 
which supports its usefulness as an indicator of the collision geometry. A max- 
imum impact parameter 6 max = 12 fm ±10% was deduced from the measured 
integrated beam intensity and the target thickness. It corresponds to the cho- 
sen trigger condition of at least 5 charged particles detected and, within errors, 
remains approximately constant over the covered range of bombarding ener- 
gies. With this information, and assuming the monotonic relation [19] between 
the impact parameter and E^ 2 , the data were sorted into six impact-parameter 
bins, each 2 fm wide and, altogether, covering the full range 0-12 fm. Particles 
of interest with Z < 2 were excluded from E 1 ^ to reduce autocorrelations, a 
procedure found useful for peripheral collisions. It was further requested that 
at least 45% of the total charge of the system are detected, a condition used to 
reject peripheral events in which the projectile residue had escaped detection. 

The kinetic-energy tensor, constructed from all identified charged particles, 
can be regarded as a global observable characterizing the preferential direc- 
tions of the particle and fragment emissions. The flow angle 0fl ow is defined as 
the angle between its largest eigenvector and the beam axis. Its distributions, 
analyzed with the weight 1 / sin6fl ow [20] , are found to be mainly a function of 
centrality, with pronounced peaks at finite angles appearing at smaller impact 
parameters. The mean value increases from between 3° and 6° for peripheral 
to about 30° for central collisions. The excitation functions are fairly flat ex- 
cept for the most central collisions. For b < 2 fm, e.g., the weighted mean flow 
angle increases from about 15° to 40° for bombarding energies between 40 and 
150 MeV per nucleon. 

The squeeze angle \l/ sq , defined as the angle between the middle eigenvec- 
tor of the kinetic-energy tensor and the reaction plane [21], characterizes the 
preferential azimuthal directions of emission. The squeeze-angle distributions 
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exhibit a clear trend as a function of incident energy and centrality (Fig. 1). 
The minima at tt/2, observed at lower energies and more peripheral impact 
parameters, indicate predominantly in-plane emissions. Peaks at ir/2, most 
strongly pronounced in the more central bins at the higher incident energies, 
correspond to a preference for azimuthal emissions perpendicular to the reac- 
tion plane, the so-called squeeze-out [21]. 
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Fig. 1. Distributions of the squeeze angle ^ sq for incident energies from 40 to 
150 MeV per nucleon (from top to bottom) and sorted into 2-fm-wide bins of the 
deduced impact parameter as indicated above each column of panels. 

The curvature of the distributions has been analyzed using their standard 
deviation cr(^ sq ). The expression a(^ sq ) / (n/ ^/l2) — 1 is positive for concave 
distributions of \I/ sq , negative for convex distributions, and zero for flat ones. 
The transition energies, E^ n , identified by a change-of-sign of this variable, 
have been determined by a linear interpolation between bombarding energies. 
Linear extrapolations were used to obtain estimates for the very peripheral 
impact parameter bins. 

The energies of the transition from predominantly in-plane to out-of-plane 
emissions are a strong function of the impact parameter (Fig. 2). They extend 
from 65 MeV per nucleon for central collisions up to about 200 MeV for 
the most peripheral collisions and exhibit an increasingly rapid rise toward 
the more peripheral impact parameters. Results obtained by representing the 
event with the kinetic-energy tensor or with a momentum tensor, i.e. without 
the weight factor l/2m, are identical at central impact parameters but diverge 
slightly in the peripheral bins. 
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Fig. 2. In-plane to out-of-plane transition energies, £ t ^ m , determined from the cur- 
vatures of the squeeze-angle distributions, as a function of the impact parameter 
b. The full and empty symbols indicate interpolated and extrapolated values, re- 
spectively. The full (dashed) line connects results obtained from analyses with a 
kinetic-energy (momentum) tensor for the event description. The horizontal error 
bars represent the widths of the impact-parameter bins and the vertical errors result 
from the systematic uncertainty of 6 max - 

The diagonalized tensors provide a global representation of the collective event 
properties. A much more detailed quantitative information can be obtained 
from the Fourier analysis of azimuthal distributions of the reaction products 
measured with respect to the reconstructed reaction plane and as a function of 
particle type, rapidity y and possibly the transverse momentum px- The first 
two coefficients, Vi and t>2, of the Fourier expansion characterize the directed 
and elliptic flow, respectively [22,23,24,25]. The reaction plane has been de- 
termined by several methods, including the flow-tensor method [26], the flow 
Q- vector method [27], and the azimuthal-correlation method [28]. 

The rapidity dependence of the directed-flow parameter v\ for Z = 2 particles, 
integrated over transverse momentum, is shown in Fig. 3. The present INDRA 
data are combined with the FOPI data [9], both measured for mid-central 
collisions with impact parameters of 2-5.3 fm and shown without corrections 
for the reaction plane dispersion. In the case of the INDRA data, the reaction 
plane has been reconstructed using the Q- vector method with the weights uo = 
p z , excluding the particle of interest ("1 plane per particle") and correcting 
for the effects of momentum conservation [29]. 
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Fig. 3. Transverse-flow parameter v\ for Z = 2 particles integrated over as a 
function of scaled center-of-mass rapidity for mid-central collisions (2-5.3 fm). The 
open and filled symbols represent the FOPI [9] and the present data, respectively. 
The numbers in the legend indicate the beam energies per nucleon in MeV. The 
vertical error bars result from the systematic uncertainty of 6 max . 

The slope of v\ as a function of rapidity is seen to rise monotonically with 
energy over the full range of 40 to 400 MeV per nucleon which is covered by 
the two experiments. Good agreement is observed in the overlap region which, 
e.g., may be verified at the incident energy of 150 MeV per nucleon which 
was used in both experiments. The coefficient v\ and its slope as a function of 
rapidity are practically zero at 60 MeV per nucleon and become even negative 
at 40 MeV per nucleon. This intriguing observation of a negative flow has 
already been reported for the lighter systems 40 Ar + 58 Ni, 58 Ni + 58 Ni, and 
129 Xe + nat Sn, provided the 1-plane-per-particle method was used [30]. For 
these systems, a balance energy, -Ebai, has been determined by associating it 
with the minima of the approximately parabolic excitation functions of the 
flow parameter which, in the cases of 40 Ar + 58 Ni and 58 Ni + 58 Ni, appeared at 
negative flow values. If the same parabolic scenario is adopted for the present 
case of 197 Au + 197 Au, the balance energy should fall below the observed zero- 
crossing of the slope dvi/dy for which the value 54 ± 4 MeV per nucleon has 
been obtained from interpolation. 



Values for the balance energy in 197 Au + 197 Au have previously been deter- 
mined by extrapolating from higher energies [31,32,33], and also by searching 
for the minimum of flow [8]. The extrapolations yielded values between 47 and 
56 MeV per nucleon with a moderate precision but, nevertheless, consistent 
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with the zero-crossing at 54 ±4 MeV per nucleon observed here. The excitation 
function of flow reported in Ref . [8] has a minimum at E^ai = 42 ± 4 MeV per 
nucleon but the measured slopes were exclusively positive for products with 
Z > 1, which is contrary to the present data. 

In order to identify the possible sources of the apparent disagreement be- 
tween these two measurements, the present data has been subjected to simi- 
lar thresholds and selection criteria as used in Ref. [8] and analyzed with the 
same azimuthal-correlation method of reaction-plane reconstruction [28], ex- 
cluding the particle of interest. The flow is represented by the slopes dv\/dy, 
determined by linear fits within the range of the scaled center-of-mass rapid- 
ity —0.5 < y C m/ycm S < 0-5- The results are summarized in Fig. 4, including 
data for 15 MeV per nucleon from a small data sample primarily collected for 
calibration purposes. 
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Fig. 4. Excitation functions of the mid-rapidity slopes of the v\ parameter for Z = 2 
particles and impact parameters 0-4 fm. The symbols correspond to the indicated 
methods and conditions used to extract the flow values. The lines represent parabolic 
fits. The statistical errors, only shown for 15 MeV per nucleon, are smaller than the 
symbol size at other energies. 
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With the conditions Z < 12 or Z < 6 and by restricting the data to positive 
center-of-mass rapidities, the parabolic excitation function of the flow param- 
eter reported in Ref. [8] is qualitatively reproduced (open symbols in Fig. 4). 
The latter condition has been identified as causing the sizable negative offsets 
of v\ at mid-rapidity in this data. Without these restrictions, the minimum 
disappears, and the flow continues to decrease to more negative values, ap- 
parently down to bombarding energies as low as 15 MeV per nucleon. The 
Q- vector method of the reaction plane reconstruction with the weights u = p z 
yields larger absolute flow values. The larger weights assigned to the faster 
and heavier fragments may cause a better definition of the reaction plane [34]. 

The origin of the negative directed flow observed in the present data and the 
strong effects of the acceptance and selection criteria are illustrated in Fig. 5. 
The left row of panels shows contour plots of the in-plane transverse velocity 
versus the center-of-mass rapidity. For peripheral collisions at high incident 
energy (top panels) the deflections of the projectile and target, as represented 
by the three-dimensional Q-vector, are small. The Coulomb repulsion from 
the heavy residues leads to the apparent depression of helium yields near 
the entrance-channel rapidities of ±0.28 and to maximum intensities at lower 
absolute rapidity, as discussed in Ref. [35]. The stronger deflection of mid- 
rapidity particles is evident from the rapid rise of v\ with y cm which, at y cm ~ 
0.2, starts to be modified by the effect of the projectile and target spectators. 

At 40 MeV per nucleon, and the same peripheral impact-parameters (middle 
panels of Fig. 5), the structure of V\ as a function of y cm is qualitatively the 
same as at 150 MeV per nucleon but compressed into a smaller range of ab- 
solute rapidity. In central collisions, at this energy, the distributions are even 
more compact (Fig. 5, bottom left). The flow vector indicates the mean trans- 
verse deflection of the forward-emitted part of the event. This concentration of 
mass and charge, apparently, causes the Z = 2 particles to be preferentially de- 
flected to the opposite side. The resulting flow, evaluated around mid-rapidity, 
is negative. Clearly, Q- vectors that do not include the heavier fragments caus- 
ing this deflection or shadowing will have different directions, and the flow 
measured relative to them may appear positive. 

The absolute sign of the deflection angle cannot be determined with the 
present methods. There is little doubt, however, that the positive flow above 
60 MeV per nucleon is associated with an overall positive deflection 
suit of the dominant nucleon-nucleon dynamics. Near and above the Coulomb 
barrier, positive deflection angles follow from the dominance of the Coulomb 
forces. The present observation of a negative flow below 60 MeV per nucleon 
may thus indeed indicate a transition toward predominantly negative-angle 
emissions for light products, as concluded in Ref. [8], but also consistent with 
polarization measurements for light fragments, even though for other reactions 
[36] . Thus, the commonly invoked picture of the attractive mean field globally 
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Fig. 5. Contour plots (in linear scale) of the in-plane component of the transverse 
velocity (p x /m) versus the center-of-mass rapidity y cm (left column) and the mean 
value V\ as a function of y C m 

(right column) for Z = 2 particles and three selected 
cases of incident energy and impact parameter as indicated. The arrows represent 
the directions of the three-dimensional Q-vector. 

balancing the repulsive effect of the collisions seems be too simple here to 
explain the observed sign change of flow. The main difference, as compared 
to lighter reaction systems, is the enlarged Coulomb field which not only has 
a strong impact on the entrance and exit channel trajectories [37] but also 
manifests itself in large repulsion effects (larger so-called Coulomb rings). 

The elliptic-flow parameter V2 permits a detailed study of azimuthal emissions 
including their mass and pp dependence, as reported for incident energies of 
90 MeV per nucleon and higher in Ref. [15]. The present data set, due to a 
higher efficiency for fragments, also reveals a significant effect of the fragment 
charge. Over the range of atomic numbers Z < 2 to Z = 9, the transition 
energy E tran , evaluated in the laboratory frame, decreases from about 100 to 
65 MeV per nucleon. The restriction to light particles, on the other hand, 
permits a comparison with other data sets and an extension of the measured 
excitation function deep into the relativistic regime. 

The excitation function of the elliptic flow for Z < 2 particles, in the rotated 
reference frame and for semi-central collisions of 4-6 fm, as measured by the 
Plastic Ball [21]), the FOPI [15], and the INDRA collaborations, is shown in 
Fig. 6. For this purpose, the present data has been analyzed with the same 
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Fig. 6. Elliptic-flow parameter V2 at mid-rapidity for Z < 2 particles from 
mid-central collisions, in the rotated reference frame. The dots, circles, and crosses 
represent the INDRA, the FOPI [15], and the Plastic Ball [21] data, respectively. 
The errors shown for the INDRA data are mainly systematic and caused by the 
uncertainty of 6 max - 

method as the Plastic Ball data, i.e. by using the kinetic-energy tensor for 
the reconstruction of the reaction plane and excluding the particle of interest. 
None of the three data sets were corrected for the reaction-plane dispersion. 
This correction is not expected to substantially alter the value of E tran (cf. 
Ref. [15]) but will affect the values of v-i at the lower energies at which the 
reaction plane is less well defined because of smaller particle and fragment 
multiplicities. 

The observed agreement between the different experiments is very satisfac- 
tory. In particular, the value of the transition energy of around 100 MeV per 
nucleon, obtained by the FOPI Collaboration, is confirmed by the present 
data. Because of the Z dependence, it is approximately 20 MeV higher than 
the global transition energy derived from the kinetic-energy tensor for the 
same impact-parameter bin (Fig. 2). Together, the three data sets constitute 
a coherent systematics of v 2 which can be expected to serve as a valuable 
constraint for transport models. 

In summary, the presented analysis of 197 Au + 197 Au collisions at intermedi- 
ate energies extends the existing excitation functions of the flow parameters 
V\ and t>2 from the relativistic regime down to the Fermi-energy regime. The 
observed transition from predominantly in-plane to out-of-plane emissions at 
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about 100 MeV per nucleoli for Z < 2 particles has been confirmed. The di- 
rected transverse flow changes its sign below 60 MeV per nucleon and becomes 
increasingly negative at lower bombarding energies. Negative flow and the ab- 
sence of a parabolic excitation function of directed flow in 197 Au + 197 Au were 
shown to be connected to the large Coulomb effects in this reaction system. 
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